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I. INTRODUCTION 



We have studied the influence of a particular kind of phason-defect on the Landauer resistance 
of a Fibonacci chain. Depending on parameters, we sometimes find the resistance to decrease upon 
introduction of defect or temperature, a behavior that also appears in real quasicrystalline materials. 
We demonstrate essential differences between a standard tight-binding model and a full continuous 
" " " , model. In the continuous case, we study the conductance in relation to the underlying chaotic map 

■ and its invariant. Close to conducting points, where the invariant vanishes, and in the majority of 
| cases studied, the resistance is found to decrease upon introduction of a defect. Subtle interference 

effects between a sudden phason-change in the structure and the phase of the wavefunction are 
also found, and these give rise to resistive behaviors that produce exceedingly simple and regular 
patterns. 

According to recent experimentsB, quasicrystals have curious (for metallic materials) transport properties. For 
example, anomalously high values of the low-temperature resistivity have been reported, and the resistivity goes 
down with introduction of defects or with increase of temperature. The anomalously high resistivity has been partly 

■ accounted for by the existence of a pseudogap at the Fermi energy. For the temperature- or defect-dependence there 
' is no clear consensus yet. In numerical simulations of 2-d systems, fluctuating behavior for the resistance of a defected 

ON (withr-fandom disorder) with respect to that of a pure system is observed, namely the resistivity can either go up or 
downa. 

Here we attempt to address these behaviors for a 1-d system by focusing on the electronic phase relations in real 
space that actually determine the conductance. We study, for example, the scattering of electrons on a quasiperi- 
odic (Fibonacci) arrangement of S- function potentials (a continuous model with the full phase coherence included). 
After determining the Landauer conductance of a finite part of such a system, then we introduce a particular 
type of defect of a step-form in hyperspaceo (called "phason-defect" from now on) and we study how this defect 
influences those phase relations. This is therefore a study of the effect of an abrupt "phase" -change in the structure 
(with "phase" defined in hyperspace) on the phase-coherence of the electronic wavefunction in real space (with phase 
• i-h , defined in the usual quantum-mechanical way) . (The role of the initial "phase" of pure (undefected) chains on the 
' resistance has been studied by other authorso). _ 

The quasicrystalline structure can induce exotic (critical) statesQ that are intermediate between localized and 
extended states. Especially in the pure Fibonacci chain all states are critical independent of values of parameters. 
One consequence that we observe, for a finite part of the Fibonacci chain, is a rather unpredictably fluctuating 
variation of the Landauer resistance with the length of the sampleu (except when we are on special energy-regions 
of intcgrability, as we will see below). But in addition to this irregular variation for a pure Fibonacci chain, we here 
also attempt the first study of the variation of the Landauer resistance for the phason-defected system as well, both 
with respect to the length of the sample and also with respect to the position of the defect. The comparison of the 
resistances between the pure and the defected system shows some interesting regular patterns, provided lhat we are 
close to special points associated with the vanishing of the invariant / of the underlying dynamical mapD (on which 
points, both problems are integrable). In addition, in the majority of values of the parameters, the pure system is 
found to be more resistive than the defected one. 

It is also demonstrated that our use of a continuous model is critical in obtaining our results. As we show, these 
results would have actually been missed had one first made a tight-binding approximation (i.e. keeping only nearest- 
neighbor overlaps). Our work is therefore a concrete example of the dangeru that the usual discretized approximations 
incorporate, especially when systems with quasiperiodicity are considered, where subtle interference effects between 
the phase of the wavefunctions and the "phase" in hyperspace can be expected. 
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II. TIGHT-BINDING MODEL 



Let us first use a tight-binding model, in order to compare with the fuller treatment that is given later. For simplicity 
we fix the Fermi (or incidence) energy in the center of the band of the outside conductor (coupled to the right and left 
of our finite system). We then study the scattering of an electron coming from the left with a wavevector k, and we 
determine the Landauer resistance^ of our finite system (ratio of the reflection to the transmission coefficient). The 
method used for this is the standard transfer matrix methodEl. 

The Hamiltonian for a general (mixed) tight-binding model is 

H = y^e n | n >< n \ + t n (\ n >< n+ 1 | + | n+ 1 >< n \ ] 

n 

with — to = t out and E — e out = 2t ou t cos ka, where k is the incident wavevector and a is the lattice constant of 
the outside conductor ("out"). The Schrodinger equation then reads in matrix language 

( * n+2 *„+l ) = M n . ( ) 



M„ = (^ -£l 0). 

The Landauer resistance (for spin-^ fermions) is then, in atomic unis (nh/e 2 ) 

Y = I T% f (1) 

with the matrix T/v = -1 .S' -1 .P/v.iS', where P is the total transfer matrix, namely P/v = -M/v__i..M/v_2...-Mo, and 
6 and S are defined by 

9 = ( e~ lkNa 00 e lkNa ) 
S=(e~ lka e lka l 1) 



A. Example : Periodic system 

In this case t\ = ... = t^-i = = t and e± = €2 = ... = ejy = e. Then P/v = M , and two cases appear naturally : 

1 sin ka 

— — —with 4> defined by (3) 
sin </> 

, E-e 



(5) 



2 

E-e , K V 2t J sinh 2 N6 

> 1 > = — r 

2£ T sin ka sinh 



b) I — r— \>l=>-= = 5-^ 5— - (6) 



with <f> defined by cosh^ = — — — (7) 

(8) 
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Case (a) gives oscillatory behavior of the resistance with size, while case (b) gives exponential increase for large 

size. In the intermediate case fc^ — ► 1 we have ^ = sil ^ fca N 2 , i.e. the resistance grows as the square of length. 
In the limit TV — > oo case (a) corresponds to allowed energies (bands), while case (b) to forbidden ones (gaps). 

One can also trivially solve a periodic system with more complex unit cell (e.g quasicrystalline approximants) . In 
this case the phase <f> depends on the internal structure of the unit cell and is a more complicated function of various 
parameters, but straightforward to determine for a given unit cell. 

Finally, the introduction of even a single defect leads to interference effects that can easily be studied through this 
method. It is foundll3 that the behaviour in the presence of a defect seems to be correlated with whether we are in 
case (a) or (b) above. 



B. Fibonacci and a Phason-defected system 



Take for simplicity an off-diagonal model (e„ = e ou t — 0). Also set E = 0. Hopping elements t's are forced to take 
two values t a ,tb arranged in : 

1) a Fibonacci way 2) with a phason-defect (corresponding to a "phase" -change of 2/t in hyperspace , with r = v/ ^ +1 
the golden mean), but in both cases in such a way as to comply with boundary conditions = to = t out = t^. This 
imposes restrictions on the possible values of N : N=5, 13,18, 26, 34, 39, 47, 52, 60,... (Their difference is seen to be a 

Fibonacci arrangement of the numbers 8 and 5). For example, for N — 39 we have 

Pure : (b)abaababaabaababaababaabaababaabaababaa(b) (9) 
2/r — shifted : (b)aababaabaababaababaabaababaabaababaaba(b) (10) 
Defected : (b)abaababaabaabahhabaabaababaabaababaaba(b) (11) 

(12) 

where in the last line we show only one of the possible places where the defect (i.e. an abrupt "phase" -change of 2/r 
in the structure; see Fig.l) can be introduced. 

We now give an analytic solution for the Landauer resistances of the two systems (pure and defected) : 
For a given N we find that 

(I) = (T- 1 -->-»)■ (14) 

\T J defected V 2 / V ' 

(15) 



with 7 = I st and s an integer power (of both signs) that is an irregularly fluctuating function of the length N (around 
the value s = corresponding to complete transparency). We conclude (for 7 > 1 ) that for s > — » I ^ ) < 

V / defected 




. We see that the appearance of such cases (which we will call "favorable" ) depends on the behavior of the 

pure 



pure chain with length (this behavior fluctuates in an irregular fashion). 

For big_szalues of N we get 50% of this favorable behaviour. (We have also carried out cases of multi-defected 
sequencestj, where we also find irregularly fluctuating behavior, but with the scale of fluctuations being much higher). 



III. CONTINUOUS SCHRODINGER EQUATION MODELS 



In this continuous case, which has the advantage of not suffering from truncation approximations, the matrix form 
of the Schrodinger equation results from matching wavefunctions and derivatives at scattering points. This yields a 
generalized Poincare mapE3 of the following form 
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with 



Mr, 



In the special case of a scattering potential of the form V(x) = ^ n A n 5(x — x n ) and with fc = w one obtainsl 



A 2 777 

2fii(ra + 1) = cosfc(x„ + i - x„) + " 2 sinfc(ir„ + i - sc„) 



r , , , ,> sinfc(a;„ + i - x„) 
A 12 (n + 1) = 

2 777 n 

^22(« + 1) = cosfc(x„+i - x n ) + ——— — T sinfc(x„+i - x„). 

2k h 



A. Comparison with Section 2 



Let us pause for a moment to see how a tight-binding approximation is usually derived: To change to a localized de- 
scription one typically writes the wavefunction ^(x) in terms of localized states 4> n , namely ^{x) — C n <j) n {x — x n ) 
with 4> n {x — x n ) = \f\^e~ Xn \ x ~ Xn \ and then neglects the overlap between distant 0„'s; assuming that only nearest- 
neighbor overlaps are significant one obtains a tight-binding model with site and hopping elementsEZl 



.1 



It is important to note then that in the special case of a continuous model with identical strengths A and quasiperiodic 
arrangements (Ax's taking two values in a Fibonacci sequence), which we later analyze in detail (see Subsection B 
below), we would just obtain a simple off-diagonal tight-binding approximation (since in this case all the e„ are 
identical). This demonstrates a defficiency of the tight-binding approximation, on which we now elaborate: 

In the case of Fibonacci arrangements (let us say of both e's and t's) and for N — F n (a Fibonacci num- 
ber) we have in the tight- binding approximation (with definitions P n = Mp n ...Mi, with M's given in the well 
knownla recursive scheme P n+ \ = P n -i.P n , with starting matrices P\ = ( E ^ a —11 0) , and P 2 = 
( 1^ — f^-l 0) . ( E ^ b — ffl 0) and hence the usual trace mapO x n +i = 2x n x n -i — x„_ 2 (with the definition 
x n = \t r {Pn))- This map has the well-known invariant / = x n+1 + x n + x 2 n _ 1 — 2x n +\x n x n -\ — 1. Straightforward 
evaluation for the above case yields 

/=^(e -e 6 )(^p- 

In our case of a simple off-diagonal model (i.e. e a = e^) (but also for the case of diagonal models (i.e. t a = tb)) 
I is E-independent as seen from , always positive, and it never vanishes. This is a deficiency of the approximation, 
as discussed below. (Note that even in the case of a mixed model, where I vanishes only for a single value of E, the 
elementary matrices M do not commute at this value, which is a major difference with the full continuous model, as 
will be discussed in Subsection B). 

We will see below that, because of implicit truncation errors, the tight-binding models miss interesting patterns 
associated with the zeros of the invariant I of the underlying dynamical map, where the basic matrices commute. 
These are actually relevant to conduction and are treated next (in Subsection B) in the continuous Schrodinger 
formulation. 
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B. Example : Periodic System 



Let us first discuss a problem where the scattering potential is an array of equally-spaced (with length a) ^-functions 
of equal strengths (A). Then, again two cases appear naturally (we take ?pr = 1, or, equivalently, replace everywhere 

^ by k * Jj- in what follows) : 

a) For | cosfca + ^ sinfca |< 1 = (^r) 2 x S ' s " n a^f with cp defined by cos0 = cos ka+^ sinfca (which is seen 

to be the same as in the usual treatment of the Kronig-Penney modeJl! for the case of allowed bands (for N — > oo)). 
We note again the oscillatory behavior of the resistance with length. 

It is easy to show that the naturally appearing phase <j) determines the crystal momentum q through (f> = qa (with 
such an identification, the resulting wavefunctions indeed satisfy the Bloch theorem, as can be easily checked). We 

can also see that the zeros of HjiJ determine the allowed values of q (from band theory) for the infinite periodic 

system. This is an important point, because it motivates our later treatment of a Fibonacci problem, where we will 
enforce the vanishing of the resistance in order to simulate the thermodynamic limit. 

b) For | cos ka + ^ sinfca |> 1 => (jf^ = (^f x ^"^/^ with (j> defined by cosh0 = cosfca + ^- sinfca. This 

corresponds to the case of gaps in the limit N — > oo. Indeed the resistance grows exponentially as N — > oo. In this 
case, wavefunctions decrease with a characteristic length A = j^^pr^ Once again the study of complex unit cells is 

of course possible. Also the introduction of a single defect results in interference effects that can be easily studied 
through this method. 



C. Schrodinger equ. with Fibonacci and Phason-defected 5-potentials 



In the pure case a recursive procedure of the form P n +i = P ra _i.P„ can also be established for this case (with 
starting matrices depending on the description). For example, in the case of equally spaced (a) (^-potentials with 
A n = {A a ,Ab} in Fibonacci arrangement, a natural description is the one given in the previous section (in terms of 
matrices M n ) (with P„ = Mp n ...Mi). Then the invariant is I = ( Aa ~ Ai> | fc | m ( fca ) [ e I > 0. (Its zeros constitute a 
periodic pattern k s = — , (n ^ 0)). 

Our system: The system, however, we will analyze in detail corresponds to (^potentials with equal strengths (A) 
but quasiperiodic arrangements {(x n — x n _i)} = {a, 6} in Fibonacci arrangementH^liJ, as already mentioned earlier. 

In this case, an alternative description (in terms of the coefficients A n and_B„ of the two linearly independent plane 
waves in the region between two <5 potentials) is more natural. One obtainsEHI 



with 



(A n+1 B n+1 ) =A(n).(AnB n ) 



A(n) = ( (1 — ih.^ e ik i x n + l- x -n) _i\_ e ik(x n+1 -x n ) i\_ e -ik(x, l+ i-x n ) ^ _|_ iA ^ e -ik(x n +l—x„) ^ 



for (^-potentials (once again we have set 2m/ h 2 = 1). In this description the Landauer resistance is given again by 
Y = \Pi2\ 2 with P the product of A's, namely P n = A(F„)...A(1). In this case the invariant / of the underlying 

dynamical map turns out to belli / = A s '" 4fc ^ a ~ b ' > i-e. also / > 0. Its zeros now constitute a periodic pattern k s = 
and will be the focus of our work in what follows. 

These special points k s are missed in the corresponding tight- binding model, as we showed earlier. In the present 
continuous model they turn out toucorrespond to commuting consecutive matrices, i.e. [P n+ i, P n ] = (which is 
consistent with the known relation^ 4/ + 2 = tr(P n .P n+ i.P~ 1 .P~_2 1 )), but even stronger, they lead to commuting 
elementary matrices ??. Indeed, it turns out that 

[A (As = a), A(Ax = &)] = A e '^° 2 - (1 - e 2 ' fc(fe ~ a) ) ( A A — 2ik — A - 2ik -A) 
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which vanishes for all special points k — {k s }. 



One can therefore say that, from a conduction point of view, the problem looks formally similar to a periodic 
problem. This is more concretely described through the following properties of the special points {k s }: 

1) They are the conducting points (extended states) surviving^ in the limit N — > oo. 

2) The Landauer resistance can be written exactly in closed form for k — {k s }, and the expression looks like that 
of a periodic system (see below). „ 

3) The set of points k s is robustQ against disorder (essentialy because of the above mentioned commutations). In 
fact it exists either for a periodic or for a random system (of two letters) . 

In what follows, we focus on these points k s — with a = r = , 6=1 (Fibonacci chain): Exactly on 
those points we get the Landauer resistance in closed form (for the Fibonacci or any disordered system of {a, b}): 

K _ ( A \ 2 sin 2 iV0 
T lfc =^ ~~ \2k~s) sm 2 <j) 

with 4> defined by | cos0| = | cosk s a + sin/c s a| = | cosfc s 6 + ^- sink s b\ (for 2m/h 2 = 1). 

Motivated by our earlier discussion on a finite piece of a periodic system (where we saw that the values of k that 
made the oscillating resistance vanish correspond to the allowed states of the infinite system), we now enforce a 
conducting behavior (resonance): For any fixed N, we choose A in such a way as to have a vanishing resistance 
(exactly at k = k s ), namely y| fc _ fc = (for both the pure Fibonacci and the phason-defected system). (This gives 
the representative behavior in the thermodynamic limit where this vanishing is actually expected for k = k s and for 
any A; see below). As a consequence we get (N — 1) different appropriate values of A (both positive and negative) 
given by 

2fc s (cos (f) s — cos k s b) 
sin k s b 

with the internal phase <p s defined by <f) s — ^ with m = 1, JV — 1 that correspond to vanishing Y\ k - k • ^ e 
(N — 1) values of the phase <j) s are symmetrically placed around 5 and as N increases they cover densely the entire 
upper half of the trigonometric circle. Consequently, the vanishing of y\ _ is representative of the behavior in the 
thermodynamic limit for any arbitrary value of A, and that is the reason we enforce it even for a finite system. By 
doing so, we then study the behavior around the conduction points k s = -^j. Below we summarize our results on 
the behavior of the resistance, of both the pure and the defected system, as k varies in the local neighborhood of the 
conduction points k s : 

1). Smallest system (consistent with boundary conditions): N = 5 

The results are given in Fig. 2, where we show the local regions around the lowest conducting point (n = 1) for one 
value of cj> s , showing that the resistance of the defected system always rises lower than the one of the pure system 
(this is also true for any <p s ). We also show the complicated behavior of the Landauer resistance in more extended 
regions around a special point and for one particular <fi s . 

Of course the problem is analytically solvable. By way of an example we give the analytical solution for n = 1 
(k s = ^rzj), m = 1 {4> s = f ) which is 

= 45.3087 (k-^j) +-.. 



-I 



30.7868 



t - 1 



This is the typical behavior for all points examined (1 < n < 5, all m's) i.e. the Landauer resistance of the 
phason-defected system around the conducting points k s , rises (from zero) always lower than the corresponding one 
for the pure system. 
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2) .iV = 13 

Extensive numerical results have been obtained in this case, where we have 3 possible defected systems, for 1 < n < 5 
and for all 12 phases <fr s . Out of those 12 cases, in 8 of them all three defected systems have Landauer resistances 
lower than that of the pure system. In only 2 cases one defected sequence rises slightly higher and also in 2 other 
cases another sequence also rises slightly higher than the pure system. We call these few cases (where the defected 
system is more resistive) "unfavorable cases" . 

3) Long chains 

We have carried out numerical calculations on chains up to more than 3000 sites, taking as a numerical convergence 
criterium the unitarity condition of our total transfer matrices. 

Because of the full phase coherence, the transmission behavior of the defected chains depends on the discrete position 
ip of the defect in interesting ways. In long chains the variable xp becomes quasi-continuous, and the discrete and 
asymmetric patterns found for small chains (see for example Fig. 3 and Fig. 4-c), now appear smooth and symmetric 
(Fig.4-a and Fig.4-b). Moreover, their actual form depends on the value of <j> a> which can be considered as a discrete 
label parametrizing the family of potential strengths A that enforce vanishing of resistance on specials points. The 
results show interesting patterns summarized below and shown pictorially in Fig. 5. 

We observe a symmetry in the qualitative behavior around the v alue <p s — \ (although the values of A corresponding 



to symmetric values of <f> s are always different, as seen from equ. Ill C . Up to this symmetry, we also observe a type 
of recurrent simplicity in the resistive behavior of the defected system, which we loosely call "cyclic" (with period 
P). For small to medium values of <j) s (in more detail, for m < N — 4P, with m defined by <f> s = TUj? and with P 
being the number of possible defect-points (note that for a given chain-length N there is a unique P always satisfying 
the inequality N < 5P for any N > 5), we observe exceedingly simple oscillatory patterns (Fig. 4). Furthermore, we 
observe that in these cases the defected system is always less resistive than the pure one. 

A qualitatively similar behavior is observed in the mirror symmetric region (see the two shaded regions in Fig. 5). 
It is interesting to point out that the above "favorable" behavior is valid for sufficiently strong absolute values of A, 

namely |A| > Ao with Ao = ^zj ~~ ["^ S^JT ' ~ 1 w ^ ^ e l° wer s i§ n corresponding to the 1st and the upper sign 
to the 2nd quadrant of the upper half of the trigonometric circle. As m increases towards P we observe additional 
fluctuations resulting in a rather regular modulated pattern, containing a finite and small number of harmonic modes 
(Fig. 6) The above mentioned "cycles" in the recurrence of the simplest behavior, occur whenever m approaches nP 
(again up to the symmetry around <p s = 7r/2). As we cross the cycles we also observe abrupt changes of ordering in 
the resistance of the defected with respect to that of the pure chain (Fig. 7). Finally, close to the symmetry point 
<t>s = § we observe largely fluctuating and modulated patterns that look self-similar (Fig. 8). These resistive patterns 
become chaotic (Fig. 8(c)) when we move sufficiently far from the special integrability points. 

(The above numerical observations are stable upon approaching the special point (i.e. decreasing e) up to the lowest 
value of e where our convergence criterion is satisfied. 

We view the above results (summarized pictorially in Fig. 5) as interference effects between the abrupt "phase" - 
change in hyperspace (that occurs at some point in physical space), and the phase of the wavefunction at that point. 
The Landauer resistance viewed as a function of the discrete variable x p (that actually corresponds to a family of 
different physical systems) is here seen to carry a "memory" of the hyperspace by showing a type of "hyper-coherence" , 
that depends on physical parameters in rather subtle ways. 



IV. COMMENT ON INTRODUCTION OF FINITE TEMPERATURE T 



The Landauer resistance at finite temperatures is given byEl 



T f(-§t)T(E)dE 



with f(E) = 1+e „k_ M) , 0=t^t- 



For N 3> 1 the transmission coefficients T(E)'s are dominant!!! in the local regions around the special points k s 
(studied on last pages). Because of the smearing of the Fermi function at low temperatures, the behavior of the 
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Landauer resistance with temperature depends on the location of the Fermi energy /x. If fi is close to a special energy 

(i.e. if fi ~ ~2m') then the decrease or increase of ^ of the defected system compared to the pure system, depends 
on whether this particular k s is a favorable or unfavorable point in the sense discussed on the previous pages. We 
conclude therefore, from our results above, that, statistically speaking, in the majority of cases the resistance will 
decrease upon increase of temperature. 



V. CONCLUSIONS 



Motivated by the anomalous transport properties of quasicrystals with defect and temperature, we have analysed 
the role of phason-defects on the Landauer resistance of a finite Fibonacci chain. In a tight-binding model with 
incident energy being fixed at the center of the band, the resistance is modified in an irregular fashion and in 
accordance with the behavior of the pure Fibonacci chain with the size of the system, with both positive and negative 
effect appearing statistically in equal percentage. In a fuller treatment of the continuous Schrodinger equation with 
(5-function potentials, the modification of the resistance depends on the location of the Fermi level with respect to 
special energies corresponding to extended states, and also on the position of the defect. For finite chains, the majority 
of cases studied show a decrease of the Landauer resistance upon introduction of defect and temperature. 

By comparing the two models, we have noted that continuous models, being free from truncation errors, preserve 
the phase coherence important for conductance and can reveal subtle additional effects. Such effects of interference 
between this coherence and the hyperspace construction have also been presented. It is interesting that most of the 
time, these give rise to exceedingly simple behaviors. 

Although we do not study the thermodynamic limit with any rigor, the finite N results are locally representative 
because of the self-similarity in the resistance v.s. length behavior. Our focus on the local behaviors around the 
special conducting points k s is justified by the fact that these are the only points relevant to conduction for N » 1, 
so that our results on resistance patterns should be representative for big systems. Some of these results could be 
tested experimentally in real chains that can be manufactured rather easily with recent advances in microfabrication 
techniques^!. 
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VI. FIGURES CAPTIONS 



FIG.l Hyperspace construction of the standard Fibonacci chain (top) and a 2/r-shifted chain (bottom) (2/r is the 
"initial phase" 9q (in units where the width of the window is 1) and it is equal to the displacement of the physical line 
in the perpendicular direction). The resulting chains, as well as the way that the defect is introduced, are illustrated 
on the top left, through a step-modification of the physical line. 

FIG. 2 (a) Typical local behavior of the Landauer resistance of the defected chain (dashed line) compared to a pure 
chain (solid line), around a special point. The defected rises from zero lower than the pure system, (b) The same 
behavior but in a more extended region of fc-space, showing various crossovers, corresponding to chaotic behavior 
shown in Fig. 8(c). 

FIG. 3 Comparison between resistances of a pure small chain (TV = 34) and all possible defected chains of the same 
length, for a fixed k close to a special point (k = k s + e). The 8 possible defected chains are shown on the bottom. 
We note a simple behavior with the defect-position (xp) and a resistance that is always lower than that of the pure 
chain (horizontal line). This simple (but discrete and asymmetric) behavior becomes smooth and symmetric for long 
chains (see Fig. 4) where xp becomes quasi-continuous, e = 1CP 4 (in units 6=1). 

FIG. 4 Corresponding results for long chains (a) for TV — 2000 the simplest possible pattern appears for (f> s = j^, 
and it is a smooth and symmetric analog of Fig. 3. This is always the pattern appearing for any TV (and m = 1). It 
also turns out that all the defected chains are always less resistive than the pure one (horizontal line). e = 10 -9 . 
(b) Corresponding result for m = 5. Again, a simple oscillatory pattern appears (with the number of bumps equal 
to m), which is generally valid for small m's. These simple patterns are quasi-continuous analogs of asymmetric and 
discrete patterns that are found for small chains, as shown in (c) which corresponds to TV ~ 200. 

FIG. 5 Global view of the resistive behavior in <j> s -\me. See text for details (end of Section 3). Shaded areas show 
regions of values of <j> a (and correspondingly of the potential strength parameter A) where the resistive behavior is 
very simple. The following three loops show the places where one recovers simple resistive patterns and transitions 
in the resistive behavior. 

FIG. 6 Modulated resistive behaviors for m between TV — 4P and P. The defected system is, in the majority of 
cases, less resistive than the pure one (horizontal line). 

FIG. 7 Recurrence of simple behaviors as m crosses P (for TV = 2000, that yields P = 472 possible defect positions). 
Note a transition in the ordering of resistances between defected and pure chains: solid curves correspond to m = P: 
dashed to m = P + 1; dot-dashed to m = P + 2. (Flat lines always give the values of the corresponding undefected 
chains). 

FIG. 8 Modulated and self-similar resistive patterns (close to special points) for chains with length TV = 3207 (which 
yields P — 757 possible defect positions), (a) For m = 1603, i.e. for <p s in the first quadrant and just before -|. (b) 
For m = 1605, i.e. for <p s in the second quadrant and two units after ^. Both (a) and (b) are given for e = 10~ 4 . (c) 
As e increases (e = 0.5) we recover chaotic behavior (the values of the parameters are the same as in case (b)). 
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